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Abstract 

Classical mechanics is formulated in complex Hilbert space with 
the introduction of a commutative product of operators, an antisym- 
metric bracket, and a quasidensity operator. These are analogues of 
the star product, the Moyal bracket, and the Wigner function in the 
phase space formulation of quantum mechanics. Classical mechanics 
can now be viewed as a deformation of quantum mechanics. The forms 
of semiquantum approximations to classical mechanics are indicated. 

While our understanding of the relation between quantum mechanics and 
classical mechanics has steadily increased over the past 75 years, as a result of 
many studies from various points of view (see [TJ- [TJ and references therein), 
few would claim that it is complete. Meanwhile, increasing attention has fo- 
cussed on the interface between the quantum and classical domains, because 
of advances in experimental science and engineering, and the associated de- 
velopment of 'nanotechnology.' 

Classical mechanics is usually formulated in real, finite-dimensional phase 
space; quantum mechanics in complex, infinite-dimensional Hilbert space. 
However, a completely equivalent reformulation of quantum mechanics in 
phase space is known |B]-[IE], which shows that quantum mechanics is a 
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deformation of classical mechanics [U] , and which provides a natural setting 
for the formulation of semiclassical approximations fUl H] . These allow us 
to explore the interface between the two forms of mechanics when approached 
from the classical side. 

In his remarkable 1946 paper, Groenewold ^T] indicated the alternative 
possibility of reformulating classical mechanics as a quantum-like theory, with 
a quasidensity operator which is not positive-definite, although few details 
were given. This idea has since been commented upon ^3] and explored in 
different ways [IS], notably by Muga et al. [17J. It is not to be confused 
with the approach to classical mechanics in the real Hilbert space of square- 
integrable phase space functions, initiated by Koopman |18j . Here we show 
that, just as quantum mechanics can be reformulated in phase space, so clas- 
sical mechanics can be reformulated in complex Hilbert space, in such a way 
that classical mechanics is seen as a deformation of quantum mechanics. And 
now there arises the possibility of exploring the interface between quantum 
mechanics and classical mechanics from the other side, the quantum side, 
with the systematic development of semiquantum approximations. 

We limit discussion to a system with one linear degree of freedom. All for- 
mulas below can be generalized to many (possibly infinitely many!) degrees 
of freedom. Our presentation is formal and heuristic; there is no attempt at 
mathematical rigor. 

A conservative classical system is usually described in terms of functions 
(classical observables) A c (q,p) on phase space, together with a probabil- 
ity density p c (q,p,t), characterizing the state of the system at time t, with 
evolution equation 

dp 

-j^- = {He, P c }p =H c Jp c , 

j d L d R d L d R 
dq dp dp dq 

Here He is the Hamiltonian function, {A,B}p denotes the Poisson bracket, 
and the superscripts L and R indicate the directions in which the differential 
operators act. The expectation value of the classical observable Ac(q,p) at 
time t is 

(A c ) (t) = J A c (q, p)p c (q, p, t) dqdp . (2) 

In (j2J) and below, integrals are over all real values of the variables of integra- 
tion. 

A conservative quantum system is usually described in terms of a complex 
Hilbert space of square-integrable state functions ip(x). Quantum observables 
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are linear operators Aq acting on state functions as 



(A Q V>)0) = J A QK (x,y)tP(y)dy, (3) 

where Aqk(x, y) is a complex-valued function, the kernel of Aq. In particu- 
lar, the canonical coordinate and momentum operators q and p have kernels 
xS(x — y) and —ih5'(x—y), respectively, where S(x) is Dirac's 'delta function.' 
If the observable quantity is real, the corresponding operator is Hermitian : 
AQx(x,y) = AQx(y,x)*. An important example is the quantum density 
operator p Q (t), which has a kernel 

p QK ( x ^y^) = XXVvOM)Vv(2/,*)* (4) 

r 

when the system is in a state described by the 'mixture' of orthogonal state 
functions tp r (x,t) with associated probabilities p r at time t. The quantum 
density operator is positive definite, with unit trace, and the expectation 
value of the quantum observable Aq at time t is 

(A Q )(t)=Tr(A Q p Q (t)). (5) 

The evolution equation for p Q (t) is 

where Hq is the Hamiltonian operator, and [Aq,Bq] denotes the commuta- 
tor. 

In order to map the Hilbert space formulation of quantum mechanics into 
the phase space formulation, the Weyl-Wigner transform W is introduced. 
For each quantum observable Aq with kernel Aqk{x,v), a corresponding 
function Aq = W(Aq) on phase space is defined by setting 

A Q (q,p) = J A QK (q-x/2,q + x/2)e ip */ n dx. (7) 

If Aq is Hermitian, then Aq is real. The Wigner density function p Q (t) = 
W{p Q (t)) I '(2'kK) is a particular case, in terms of which the quantum expec- 
tation value (JSJ) can be rewritten as 



(^q)(*) = J r A Q( ( liP)pQ( ( l,P, t ) dqdp. 



This has the appearance of the classical average (|5jL but while the Wigner 
function is real and normalised, it is not in general nonnegative everywhere on 
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phase space, and consequently can be interpreted only as a quasiprobability 
density 

In order to describe dynamics in the phase space formulation, the cel- 
ebrated star product and star (or Moyal) bracket of quantum phase space 
functions are introduced [HI HU H2j : 

A Q *B Q = W(AqBq) , 

{A Q ,B Q l = ±W([A Q ,B Q ]) 

= ±-(Aq*B q -Bq*A q ). (9) 

Then q-kp = qp + ih/2, p-kq = qp — ih/2, q 2 *p 3 = q 2 p 3 + Zifiqp 2 — 3h 2 p, etc. 
The quantum evolution (JHJ) is now replaced by 

dp Q (t) 



dt 



{Hq, Pq I, (10) 



where H Q = W(H Q ). 

For suitably smooth Aq and Bq, in particular polynomials in q and p, it 
can be shown from that 

{A q ,B q \=A q GBq, G=Hffln[5j], (11) 

where the sine function is to interpreted by its Taylor series, and J is as in 
(0). For more general Aq, Bq, such an expansion has only an aymptotic 
meaning, so that (fTUj) leads to 



dpJt) h 2 

h 4 



+ §^H Q J 5 p Q ..., (A-0). (12) 

Equations (JHJ) and (|12jl are to be compared with their classical counter- 
parts (j2J) and ((TJ), which are 'obtained' when % — > 0. It is not our purpose here 
to discuss the subtle mathematical difficulties associated with this limiting 
process PQE]. Suffice it to say that (JHJ) and (fT2"|) form a natural starting point 
for discussions of the classical limit, and of semiclassical approximations to 
quantum mechanics as % approaches 0. 

We now stand the foregoing on its head. With each classical phase space 
function A c (q,p) we associate a linear operator Ac = W~ l (A c ). This defines 
Aq as the operator with kernel 

A CK (x,y) = ^-J A c ([x + y]/2,p)e^ x -y^ h dp. (13) 
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If Ac is real, then Ac is Hermitian. This is the usual Weyl mapping [H] 
from functions to operators, but our intention here is not to quantize, but 
to reformulate classical mechanics in complex Hilbert space. It may then be 
objected that Planck's constant is not available to us in a classical theory. 
We treat H for the moment as a parameter with dimensions of action, whose 
value is to be specified at our convenience. 

As a special case, we have the Groenewold density operator [TT] 

p c (t) = 2nhW-\ Pc (q,p,t)), (14) 

This can be seen to be bounded, with unit trace, but unlike a true quantum 
density operator, it is not always positive-definite. Just as the Wigner func- 
tion p Q (q, p, t) is only a quasiprobability density, so the Groenewold operator 
p c (t) is only a quasidensity operator. But just as quantum averages can be 
calculated using the Wigner function in the 'classical' formula (jHJ), so classical 
averages can be calculated using p a (t) in the 'quantum' formula 

(Ac)(t)=Tr(A c p c (t)), (15) 

where Ac is the operator corresponding to the classical function Ac(q,p)- 
These averages do not, of course, involve the parameter h. 

In order to describe classical dynamics in complex Hilbert space, we first 
introduce a distributive, associative and commutative 'odot' product of op- 
erators, 

A c B c = W- l {A c B c ) = B C QA C . (16) 

Then for example, q © p = p © q = (qp + pq)/2, q 2 p 3 = p 3 q 2 = 
(q 2 p 3 + 2qp 3 q + p 3 q 2 )/A, etc. More generally, {q k p 1 } {q m p n } = {g fc + m ^+ n } ) 
where {q r p s } denotes the Weyl-ordered operator [EHE! corresponding to the 
classical monomial q r p s . This follows from (jlfij) because {q r p s } = W^ 1 (q r p s ). 

Most generally, it can be seen from (|13|) that the kernels of the operators 
Ac, Be and A c B c are related by 

{Ac B c )k(x, V) = (Be © Ac) K {x, y) 
= J A CK (l3x + y- 2u}/4, [x + 3y + 2u]/4) 
xB CK {[3x + y + 2u}/4, [x + 3y - 2u]/4) du . (17) 

It is helpful to introduce the notations 

A q = dA/dq , A qp = d 2 A/dqdp , . . . 
A q = hA,p], A qp = (h 2 [q,[A,p}},... (18) 
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and to note that, because A qp = qG (AGp), etc., and 

W\AGB) = ^[A,B], (19) 

we have W _1 (A 9P ) = A qp , etc. In (fT£j) . q and p are the usual canonical 
operators, except with commutator involving the parameter h, whose value 
has not yet been fixed. 

To describe classical dynamics, we need to introduce a new bracket, equal 
except for a convenient factor to the image of the Poisson bracket under the 
inverse Weyl-Wigner transform. We set 

[Ac,B c ] e = ihW-\{A c ,B c }p) 

= ih{A Cq B Cp - A Cp B Cq ) . (20) 

Now the classical evolution equation is replaced by 

■[H c ,p c ) Q . (21) 



dp c _ 1 
~0f ilr" 



We emphasize that this reformulation of classical mechanics in terms of linear 
operators on Hilbert space, incorporating the arbitrary parameter h, and with 
key equations (fT3j) and (}2"Tj) . is entirely equivalent to the usual phase space 
formulation. We can go back and forth between the two descriptions with 
the help of the Weyl-Wigner transform W and its inverse W _1 . 

Next we make an expansion of the odot bracket, analogous to the expan- 
sion (fTTj) . Noting that 

6 = sin#(l + #76 + 7fl 4 /360 - ... , |0|<tt, (22) 



we write 
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AJB = AGB + -\-\ AJ 2 GB 
6\2 



7 (K^ 



~ I AJ A GB - ... = AGB 



360 \2 / 

+- (-J (A qq G B pp — 2A qp G B qp + A pp G B qq ) + . . . 



(23) 



and then, applying W 1 to both sides, 

[A,B] Q = [A,B] 



+ g f 2 ) "^wJ _ 2[Ay p , B qp ] + [App, Bqq]) + . . . 



(24) 
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The series ()23j) and (|24|) terminate if at least one of A and B is a polynomial 
in g and p. For more general A and B, we may expect that the series have 
well-defined meanings as asymptotic expansions when H — > 0. 
The classical evolution then takes the form 

«9t ih [ c,Pcl 

ih — — — 

- . . . , (ft -+ 0) . (25) 

If He is a polynomial in q and p, then this series terminates and the asymp- 
totic result becomes exact. 

If H c = H(q,p) = p 2 /(2m) + V(q), then (j23J) reduces to 

f ->(«,». a,] -|[H4).w 

7 ? - ft 3 

-^[^ ) W),Po_] + ---, (ft-0), (26) 

which is an analogue of Wigner's equation for the evolution of his density 
function [TO] . 

If we now identify h with Planck's constant, we see that the equations 
(0) and (jUJ) of quantum mechanics are obtained formally as h approaches 0, 
and that classical mechanics can be regarded as a deformation of quantum 
mechanics, with deformation parameter h. Most interesting is that (|TH|) . 
taken with (|23|) or (|2l)|) . may be expected to form a suitable starting point 
for semiquantum approximations to classical mechanics, analogous to semi- 
classical approximations to quantum mechanics. Successive approximations 
will be associated with successively later terminations of the series (I25j) or 
(J2HJ). Note that h may then appear in the corresponding approximations to 
classical averages calculated using (JEJ), in contrast to its non-appearance in 
exact classical averages. 

These results may seem paradoxical. We have introduced h into a re- 
formulation of classical mechanics, without affecting its predictions in any 
way, and see that as that parameter approaches 0, the equations of quantum 
mechanics emerge. Usually we say, speaking loosely, that classical mechan- 
ics is obtained from quantum mechanics as h approaches 0. Viewing things 
from the perspective provided by the above results, we argue that it is more 
appropriate to say that classical mechanics and quantum mechanics become 
asymptotically equivalent as h — > : the interface can be approached from 
either side. 
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We close with a few remarks about interesting issues arising: 

• The fundamental importance in quantum mechanics of the spectra of 
selfadjoint operators, the superposability of state functions, and the nonuni- 
tary change in the density operator following a measurement, are obscured 
in the phase space formulation. They underlie the determination of averages 
(|5|) and of initial values of Wigner functions. On the other hand, the complex 
Hilbert space formulation of classical mechanics begs the question: What are 
the relevance to classical mechanics, when formulated in this way, of operator 
spectra and the superposability of complex vectors? 

• Consider a normalized classical density at some fixed time given by 

/ a/? __ (ag 2 +/3 p2) 



fhfap) = —e~ {ar+f * } ■ ( 27 ) 

7T 

The corresponding quasidensity operator p c has kernel 

Pck ( X , y) = /V^)VV(*-,)V(4/^) (2g) 



It is easy to check that this operator is bounded, with unit trace, but it is 
not in general positive definite. If a(3 = l/(h) 2 , so that the product of the 
uncertainties of q and p equals h/2, the kernel factorizes: 



o CK = MxWW , 1>{x) = { - e-« 2 / 2 , (29) 



and the operator has the form of a true, positive-definite density opera- 
tor, corresponding to the pure coherent state ip(x). More generally, a little 
thought shows that the only positive-definite quasidensity operators are those 
corresponding to convex linear combinations of Gaussian p c {q,p), each with 
the product of the uncertainties in q and p equal to h/2. At the other extreme, 
as a — > oo and (3 — > oo, then p c (q,p) — > 6(q)8(p) and p CK (x, y) — > 25(x + y). 
This defines the starting point of a classical trajectory, as described in the 
Hilbert space formulation. 

• Consider a classical system exhibiting chaos [3] , for example the Henon- 
Heiles oscillator with 2 degrees of freedom and Hamiltonian 

H c = H(q 1 ,q 2 ,pi,p 2 ) 
= a( Pl 2 +p 2 2 ) + b( gi 2 + q 2 2 ) + cq, (3q 2 2 - gi 2 ) . (30) 

This system is described in Hilbert space by ()15j) and (the obvious general- 
ization of) (|26jh with the series terminating after the terms of order h. If 
we choose a Gaussian initial density, generalizing (|27|) . with arbitrarily small 
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uncertainties in the dynamical variables then, with the help of a computer, 
we can in principle track the average evolution of the classical system, again 
with arbitrarily small uncertainties, and even if the motion is chaotic, while 
working in the Hilbert space formalism. This is remarkable because in the 
leading 'quantum approximation,' obtained by neglecting the terms of order 
h in , the classical chaos is suppressed [TJ |3] • 

• The new bracket has the 'odot derivation property' and 'odot Jacobi 
identity' which it inherits from the Poisson bracket: 

[A, BQC] G =CQ [A, B] G +BQ [A, C] e 
[[A, B) e ,C) e + [[B, C] @ , A) @ + [[6, A) Q ,B) e = . (31) 

Poisson algebras of phase space functions, and associated groups, should 
translate into interesting odot operator structures in Hilbert space. 
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